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ABSTRACT 

A previous  investigation  by  one  of  the 
present  authors  extended  the  method  of  Edge- 
Functions  to  the  determination  of  natural 
frequencies  and  associated  mode  shapes  of  free 
vibration  of  thin  elastic  plates  with  a variety 
of  boundary  conditions.  The  present  Investiga- 
tion further  extends  this  approach  to  the  case 
of  a shallow  elastic  spherical  shell  of  n-slded 
polygonal  plan-form.  Natural  frequencies  and 
associated  mode  shapes  together  with  boundary 
residuals  (indicating  how  well  the  approximate 
solution  has  satisfied  boundary  conditions)  are 
readily  displayed  by  the  computer  program 
offered.  Results  obtained  in  the  present  work 
are  shown  to  be  in  excellent  agreement  with 
existing  values  for  the  specific  boundary 
conditions  treated.  The  present  approach 
involves  only  modest  computer  efforts  but  offers 
the  significant  feature  of  rapid  determination 
of  how  precisely  the  boundary  conditions  have 
been  satisfied  along  each  edge  of  the  shell. 

INTRODUCTION 

The  objective  of  the  present  investigation 
is  to  present  a new  generalized  approximate 
approach  to  the  problem  of  free  vibrations  of 
thin,  shallow,  elastic  shells.  The  technique 
is  applicable  to  shell  structures  having  an 
arbitrary  polygonal  contour  in  the  base  plane 
and  arbitrary  boundary  conditions  along  this 
contour.  The  specific  example  of  a shallow 
spherical  shell  is  investigated  in  detail.  The 
present  investigation  employs  the  Edge-Function 
method  recently  successfully  applied  to  deter- 
mination of  natural  frequencies  of  free 
vibration  of  initially  flat  thin  elastic  plates 
by  I.H.  Tai  and  W.A.  Nash  [1].  The  study  [1] 
and  the  present  work  permit  rapid  determination 
of  plate  and  shell  natural  frequencies  (respec- 
tively) and  associated  mode  shapes  with  the 
expenditure  of  only  very  modest  amounts  of 
computer  effort.  Further,  the  present  approach 
permits  comprehensive  evaluation  of  the  errors 
involved  in  this  approximate  technique;  a 
feature  that  is  absent  in  most  other  approaches 
for  determination  of  plate  and  shell  natural 


frequencies. 

FUNDAMENTAL  EQUATIONS 

The  basic  differential  equations  employed 
in  this  work  to  describe  the  behavior  of  a 
thin  elastic  shell  In  the  absence  of  applied 
loads  are  those  derived  In  Kraus  [2].  The 
system  of  equations  reduces  to  two  simultaneous 
differential  equations  for  the  normal  deflec- 
tion w and  the  stress  function  $ of  the  form: 

2 

DvVw  + V2*  * -oh  2"  (1) 

* ar 

jL  vVq  - v2w  * o (2) 

where  v2  is  the  Laplace  operator,  k and  k are 

a y 

curvatures  in  the  x and  y directions  respec- 
tively, V2  " ky  fjr  ♦ kx  , E represents 
Young's  modulus,  h the  shell  thickness,  v Is 
Poisson's  ratio,  p the  shell  density,  t denotes 
time,  and  D * Eh2/12(l-v2).  Normal  and  shear 
stress  resultants  as  well  as  bending  and 
twisting  moments  may  be  expressed  in  terms  of 
w and  * as  Indicated  in  [2].  The  tangential 
displacements  u,  v along  the  x and  y axes 
respectively  may  be  expressed  In  terms  of  w 
and  p from  relations  presented  In  [2],  For 
example,  the  first  of  these  relations  is: 


Let  us  consider  the  case  of  a shallow 
spherical  shell  for  which  kx  • ky  * k ■ 1/R 
where  R Is  the  radius  of  curvature  of  the  shell 
middle  surface.  If  both  w and  * exhibit  a 
harmonic  dependence  on  time,  i.e.  if 


w ■ W(x,y)  sinuit  (4) 

4 * *(x,y)  sinut  (5) 

then  Equations  (1)  and  (2)  reduce  to 

0V4W  + kV2$  * ohu2VI  (6) 

- k^M  * 0 (7) 
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where  u Is  the  natural  circular  frequency  of  free 
vibration.  It  Is  possible  to  uncouple  (6)  and  (7)  Into 


4 4 

p V * * 0 


p4  ■ g(pu?  - Ek2)  (10) 

Since  there  are  four  conditions  to  be  specified  on  each 
edge  of  the  polygonal  boundary  the  solutions  for  w and 
<p  must  contain  four  arbitrary  constants. 

For  convenience  let  us  set 

i|>  * V4W  - p4W  (11) 

Substitution  of  this  relation  Into  (8)  leads  to 

V2*  * 0 (12) 

FORMULATION  OF  SOLUTION 

The  projection  of  the  shell  onto  the  x-y  plane  Is 
considered  to  be  a two-dimensional  convex  simply-con- 
nected region  R,  the  closed  boundary  B of  R being  a 
polygon  of  J sides.  As  In  Figure  1,  a set  of  rectangu- 
lar Cartesian  coordinate  axes  Oxy  Is  chosen  and  the 


rq  + l x, 


Pi  / 


.'  \ Po 


V Pit. 


Yj  A*i 


FIGURE  1 

vertices  and  sides  of  the  polygon  B are  numbered  1 to  J. 
The  typical  vertex  Pj  has  coordinates  (Xj.Yj)  in  the  Oxy 
reference  frame  and  the  typical  side,  the  Jth,  Is 
PjPjt,,  having  length  aj  and  making  an  angle  d> j with 
the  positive  direction  of  the  Ox  axis.  Associated  with 
each  side  J Is  a set  of  edoe-axes  Xjyj,  having  origin 
at  P,  and  directed  along  the  (inward)  normal  to  that 
side. 

The  fundamental  series  solution  of  equation  (12)  In 
rectangular  coordinates,  derived  using  the  separation  of 
variables  technique.  Is 

* ■ £ {[A^lntmx  ♦ c^)  + Bmcos(mx  + exile""1*  + 

(m) 

[Cms1n(mx  + a,,)  ♦ 0mcos(mx  ♦ exile*™*)  (13) 

where  m.  A,,,.  8m,  Cffl  and  Dm  are  as  yet  undetermined 
constants.  The  phase  shifts  and  B,,,.  although 
unnecessary  here,  are  Included  for  sake  of  generality 
later.  Since  the  Laplaclan  operator  T2  Is  Invariant 
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under  transformations  which  produce  translation  and 
rotation  of  axes.  It  follows  that 

(14) 

j 8*j 

The  solution  of  this  equation  will  then  take  the  form 

* * (£j){CAWs1n(mjxj  + ^ + Bmcos<"jxj  * 

6m>le""j*J  + [C^IntmjXj  ♦ o*)  ♦ ^cosdyj  ♦ 

Bm)]e+™j*j)  (15) 

In  terms  of  the  coordinates  (*jyj)  of  a point  referred 
to  the  jth  system  of  axes.  Subscripts  and  superscripts 
j have  been  added  to  m and  to  the  constants  A,,.  Bm,  Cm, 
and  Dm  to  Indicate  association  with  the  jth  side  of 
the  boundary  B.  Clearly,  similar  expressions  arise  for 
solutions  associated  with  each  of  the  J sides  of  the 
polygon  and  all  such  solutions  may  be  superposed,  the 
governing  differential  equation  being  linear.  Thus,  a 
general  solution  may  be  put  In  the  form 

* = j-1  (mj)  {[A"S,"(,Vj  + <*>  + Bicos(™Jxj  + 

B„|)]e-mj*J  ♦ [CjJs1n(mjXj  ♦ %)  ♦ ofeosim^  ♦ 

Bm)3emJ*j)  (16) 

Analogous  to  the  solution  for  a simple  rectangular 
region,  the  constants  nij  In  (16)  are  chosen  so  that 

mj  * 2tWa.j  (17) 

for  M a non-negative  Integer,  In  order  to  facilitate 
generation  of  the  boundary  condition  equations  later. 
Obviously  the  function  contributes  to  the  determina- 
tion of  the  displacement  components  and  stress  couples. 
Thus,  In  general.  In  formulating  the  boundary 
conditions  for  the  typical,  jth,  side,  the  coefficients 
In  (16)  associated  with  the  jth  side  can  be  made 
Implicitly  dependent  on  the  relevant  boundary  condi- 
tions. In  the  general  solution  at  any  point  within 
the  region  R the  displacement  components  and  stress 
resultants  will  be  composed  of  contributions  from  each 
side  of  B,  each  contribution  being  dependent  on  the 
boundary  conditions  Imposed  on  that  particular  side. 

Thus,  Invoking  St.  Venant's  principle.  It  Is  physically 
desirable  that  such  contributions  decay  with  Increasing 
distance  from  the  corresponding  side.  So,  In  (16)  the 
coefficients  CjJ  and  0^  are  set  equal  to  zero.  Conse- 
quently * Is  chosen  to  be 

J 1 l 

* * J.1!  jtj  iMSln(™jXJ  + «m>  + BlMcos(mjxJ  * 

0m)]**mJyJ  (18) 

In  which  nij  Is  defined  In  (17). 

Each  of  the  fundamental  terms  e"™j*Js1n(mjXj  + a,,) 
and  e ™J*j  cos(mjXj  + B^,)  In  <P  now  possesses  the 
advantageous  properties  of  being  directly  associated  with 
the  j^  side  of  the  boundary  and  of  decaying  In  contri- 
bution to  the  overall  solution  with  increasing  distance 


yj  Into  the  region  from  that  side.  Such  functions, 
similar  both  In  form  and  notation  to  those  employed  by 
Tal  and  Nash  [1],  are  termed  edoe-functlons  after 
Quinlan  [3],  Comparable  expressions  for  0 may  be 
written  but  are  omitted  for  brevity.  If  one  substitutes 
(18)  Into  (11)  and  also  carries  out  a comparable  pro- 
cedure for  ♦,  one  finally  obtains 

W - E Z {[-  ^ A3  e‘mJyj  + A3  e'(mj  + PZ)VZyj  + 

j-1  H ^ 1H  3,1 

A^e'^j  ” p 3 y33s1n(mjXj  + c^,)  + 

[-  \ B3Me""jyj  + B^e'^j  + pZ)VZyj  + 

P 

2 2 M2 

B^e'h  ' p > ^cosOPjXj  + Sn,))  (19) 

and 

* " j-\  N{C*  7***  ' ^ mjyjAlM)e’mjyj  + 

>V(^pZ),/VAi,e-(fp2l'/2yJ)] 

P 

•s1n(mjXj  + °*]  + C'  ^f(B2M  ' mJyjB1M)' 

• e-mJyj  ♦ + pZ)1^yJ  - 

, iJ  n2, 1 /2 

B4Me"(  i ' P ’ yJ)]cos(mjXj  ♦ Bm)>  (20) 

The  expression  for  the  transverse  deflection  (19) 
reduces  In  the  case  of  zero  curvature,  k * 0,  to  the 
corresponding  expression  obtained  by  Tal  and  Nash  [1] 
for  the  transverse  deflection  of  a thin  flat  plate. 

Shell  boundary  conditions  are  expressed  In  terms  of 
various  combinations  of  the  parameters  W,  0,  the  In- 
plane displacements  uq  and  vq  along  and  perpendicular 
to  the  q-th  side  of  the  boundary,  the  rotation  6 about 
the  tangent  to  the  q-th  side, the  rotation  Bq  about  the 
normal  to  the  q-th  side,  the  In-plane  stress  resultant 
Nq,  the  bending  moment  Mq  about  the  q-th  side,  the 
twisting  moment  Mpt  about  that  same  side,  the  Klrchhoff 
effective  shearing  stress  resultants  Tq  and  Vq,  and  the 
In-plane  shearing  stress  resultant  Njjt<  all  associated 
with  the  q-th  edge.  Expressions  for  each  of  these 
twelve  parameters  are  derivable  In  terms  of  the 
coefficients  A^  and  B3M  but  are  omitted  for  brevity. 
However,  It  Is  convenient  to  symbolize  these  boundary 
functions  In  the  form 


where  t Is  an  Integer  ranging  from  1 to  12  for  the 
above  parameters  W,  ♦,  uq  etc.  respectively.  Further, 
let  us  Introduce  the  function  A(E^  to  denote  the  edge- 
function  contribution  to  A^  stemming  from  the  J-th 
coordinate  system  and  associated  with  the  particular 
value  14.  In  this  case  each  of  the  above  boundary 
functions  may  be  written  In  the  form 
J 1 

A,  • £ l AtEi  (21) 

1 J-l  M 


It  Is  necessary  to  determine  values  for  the  fre- 
quency as  well  as  the  coefficients  a|h  and  b|h  so  as 
to  satisfy  boundary  conditions.  Thus,  boundary 
conditions  require  that  four  of  the  functions  At 
be  specified  on  each  edge  of  the  polygonal  boundary. 
This  Is  accomplished  as  follows. 

If  Pq  Is  ary  point  (xq,0)  on  the  qth  side  of  the 
boundary,  as  In  Figure  1,  the  function  At(Pq)  Is 

specified  for  four  values  of  the  parameter  t » t1, 

2 3 4 m 

tq,  tq  and  tjj,  say.  Consequently,  equation  (21) 

Implies  that 


At(V  ■ J.*,  E 


ft1  t2  t3  t4- 

1 V q’  V V 


1,2 J 


In  order  to  obtain  Identity  equations  for  the 

coefficients  A3M  and  B3M,  It  Is  noted  that  the  series 

of  terms 

I AtE2<p;> 

(M)  t M q 

In  the  summation  (22)  does  not  contain  negative 
exponential  factors  since  yq  * 0,  so  that  each  term 
Involves  only  xq.  Since  the  negative  exponentials 
tend  to  diminish  boundary  Influences  this  series  of 
terms  Is  dominant  and  consequently  (22)  may  be  written 

MV  ' J,  S AtEH(Pi)  ■S{j1(h1A1H  + 9,B1H) 

4 


s,n<Vq  + V 


rA<w  + 9,..B?m) 


cos(m  x + Bj)  (23) 

J 

The  dash  In  the  sumnatlon  Ej , Indicates  that  the  term 
J * q Is  omitted.  The  factors  h and  g In  (23)  are 
Independent  of  xq  and  may  readily  be  obtained  from  the 
expressions  for  AjEjj  by  setting  yq  * 0.  Thus  the  left 
hand  side  of  (23)  Is  some  function  H(xq)  of  xq  and  the 
Identity  can  be  satisfied  provided  that  H(xq)  can  be 
expanded  In  the  trigonometric  form  of  the  series  on 
the  right-hand  side.  On  choosing 

mq  * 2Nir/aq,  M - 0,1,2,...  (24) 


the  right-hand  side  of  (23)  may  be  considered  the 
Fourier  series  expansion  of  H(xq)  In  the  range  [0,aq]. 

Several  methods  of  approximately  satisfying  the 
boundary  Identities  (23)  may  be  employed.  Perhaps  the 
most  obvious  Is  to  multiply  (23)  by 

sin  x dx 
aq  q q 

and  Integrate  from  0 to  aq,  then  repeat  using  the 
cosine  function.  This  leads  to  a set  of  equations 
of  the  form 

G(xq)  ■ 0 (26) 

which  Is  actually  an  Infinite  set  of  equations  In  an 
Infinite  number  of  unknowns.  It  Is  of  course  necessary 
to  truncate  these  equations  at  some  level,  say  lj,  for 

379 


each  side  j of  the  polygon. 

Because  the  truncation  of  the  harmonic  series  expan- 
sion of  G(xq)  at  N 


Lq  implies  that  the  approximation 


G(V  * K + $l*coi  if 


FNsin?^Xq]  (27) 

is  being  used,  another  way  to  approximate  G(xq)  is  to 
use  trigonometric  interpolation  to  fit  a finite  trigono- 
metric series  to  G(xq)  by  a discrete  least  squares 
method.  As  is  shown  in  reference  [4],  if  a discrete 
Fourier  series  interpolation,  based  on  discrete  least 
squares  fitting  at  2k*-l  equidistant  Internal  points 
in  the  Interval  [0,aq]  is  used  in  approximation  (27), 
then 


1 

k* 

1 

k* 


2k* 

Z 

k*0 
2k*- 1 

Z 

k*l 


wkG(Xk)cos  ^ xk 
wkG(Xk)sin  xfc 


N * 0.1,2 L„ 


N * 1,2,. ...L„ 


where 


kaq/2k* 


(28) 

(29) 

(30) 


and  wk  are  weight  factors  defined  by 

wk  * 1/2:  k * 0,  k * 2k* 

* 1:  0 < k < 2k* 

The  expressions  (28)  for  EN  and  may  be  set  to  zero 
thus  generating  the  requisite  set  of  simultaneous 
boundary  equations.  This  interpolation  method  of 
setting  up  the  simultaneous  equations  has  the  advan- 
tage, over  the  first  method  above,  of  requiring  less 
computational  time.  Furthermore,  it  is  found  in 
oractice  that  there  is  very  little  difference  in 
accuracy  between  the  results  obtained  from  either 
method. 

The  set  of  equations  arising  from  expression  (27) 
in  which  the  coefficients  of  the  expansion  must  be  set 
to  zero,  ensures  that  the  boundary  conditions  (22  ) are 
satisfied  at  all  points  of  the  qth  side,  q • 1.2.....J, 
except  possibly  at  the  vertices.  If  the  boundary 
conditions  are  not  automatically  satisfied  at  the  end 
points  of  each  side  of  the  polygon,  tliey  must  be 
inposed  at  these  points.  Thus  it  is  required  to  set  to 
zero  both  the  boundary  function  and  its  derivative  along 
the  relevant  side,  so  that 


G(»q) 


and 


3 

3x. 


G<»„> 


(31) 


(32) 


aq  of  the  Interval  and  for 


for  both  ends  xq  * 0 and  xq 

each  of  the  boundary  functions  corresponding  to  t 
tq,  tq,  tq,  and  tj.  The  effect  of  these  point 
equations,  which  are  referred  to  as  vertex  equations, 
is  not  alone  to  ensure  that  the  boundary  conditions 
are  satisfied  at  the  end  points  of  the  Interval,  but 
also  to  increase  the  rate  of  convergence  of  the  Fourier 
series  from  order  1/M  to  at  least  order  1/M3.  Such 
equations  clearly  give  rise  to  the  necessity  of  intro- 
ducing additional  unknowns  in  the  basic  solutions  (19) 


mq  ■ 2Mi./aq. 


and  (20).  Since,  in  general,  there  are  four  vertex 
equations  arising  for  each  boundary  condition  on  each 
side,  a total  of  sixteen  additional  unknowns  associated 
with  each  side  are  required.  These  are  provided  by 
using  fractional-edoe-functions.  The  latter  are  edge 
functions  of  the  types  in  expressions  (19)  and  (20) 
within  the  summation  signs.  Whereas  in  the  above 
harmonic  edge-functions 

H * V °m  ’ °»  Bm  * 0 

(33) 

for  the  fractional-edge-functions  new  values  of  these 
parameters  are  defined,  so  that 

mq  * 2V/aq;  °m  * 01  * 0 (34) 

where  V„  is  non-integral.  The  sixteen  requisite 
unknowns  can  be  provided  by  including  in  the  solution 
two  fractional -edge-functions  generated  simply  by 
choosing  two  different  values  for  V^.  They  may  be 
Incorporated  in  the  general  solution  by  adopting  the 
convention  that  they  correspond  to  M * -1  and  M « -2 
in  the  sumnations  in  (19)  and  (20.  It  may  be  noted 
that,  in  theory,  the  solution  is  independent  of  the 
choice  of  the  vertex  numbers  and  of  and  ^ 
for  the  fractional-edge-functions.  However,  in 
practice  it  is  found  that  a judicious  choice  of  values 
will  Increase  the  convergence  of  the  truncated  series 
solution. 

It  was  found  in  [1]  that  the  problem  of  free 
vibrations  of  flat  plates  could  be  successfully 
treated  by  use  of  a combination  of  edge  functions  and 
fractional  edge  functions.  That  is,  the  prescribed 
boundary  conditions  could  be  satisfied  to  the  desired 
degree  of  accuracy.  In  the  case  of  shell  vibrations 
this  is  not  the  case  and  it  becomes  necessary  to 
employ  not  only  the  two  types  of  functions  used  in 
plate  analysis,  but  in  addition  a third  type  of 
function  termed  a shell  polar  function.  This  function 
contributes  to  the  time-dependent  deflections  in  the 
interior  regions  of  the  shell  away  from  the  boundaries, 
whereas  the  edge  functions  contribute  largely  in  the 
imwediate  vicinity  of  the  shell  edges.  The  shell 
polar  functions  are  of  the  form 

* - Re[EZX  ♦ F ZZX+1] 

W ■ R e[F\ZX] 

where  z * x ♦ 1y  (origin  at  any  convenient  interior 
point),  z is  the  complex  conjugate,  E,  F,  and  F'  are 
complex  constants,  and  x is  an  Integer.  Details  of  the 
derivation  and  accompanying  properties  of  these 
special  functions  are  presented  in  the  Appendix. 

The  shell  polar  functions  are  appended  to  the 
solutions  to  replace  some  (or  possibly  all)  of  the 
fractional  edge  functions.  The  selection  of  the  ratio 
of  fractional  edge  functions  to  shell  polar  functions 
is  made  on  the  basis  of  experience  and  judgment  in  use 
of  this  technique  but  with  the  requirement  that  boundary 
conditions  be  satisfied  to  the  prescribed  degree  of 


(35) 
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'accuracy.  In  the  present  Investigation  a combination 
of  half  edge  functions  and  half  polar  functions  was 
found  best  for  satisfying  the  prescribed  boundary 
conditions. 


For  the  case  of  a clamped  edge  spherical  shell  the 
first  four  natural  frequencies  are  indicated  below,  as 
well  as  those  found  by  the  Malkina  method  [6]  again 
using  the  greatest  inscribed  circular  base. 


Application  of  the  boundary  conditions  leads  to  a 


Frequency 


system  of  homogeneous  equations 

in  the  unknowns  A3M 

Present  Method 

and  B^  of  the  form 

L=1 

L=7 

6746.24 

6746.24 

6794.72 

6788.39 

6992.20  7135.09 

ICab(“)'\,  * °*  >>  = 1.2 

(36) 

Malkina  [6] 

6746.24 

6796.89 

6998.91 

: 


\ 


where  Ab  represents  A^H  and  B^.  The  Cab  are 
coefficients  obtained  from  the  Interpolation  method 
mentioned  previously  and  are  functions  of  the  unknown 
frequency  u>.  The  frequencies  are  found  by  setting  the 
determinant  of  the  system  equal  to  zero,  viz: 

|Cab(w)|  - 0 (37) 

Values  of  in  are  determined  by  an  suitable  iteration 
method,  such  as  the  bisection  procedure. 


NUMERICAL  RESULTS 

Let  us  consider  the  free  vibrations  of  a thin, 
elastic  shallow  spherical  steel  shell  of  square  plan- 
form.  The  radius  of  curvature  is  30.0  Inches,  the 
shell  thickness  0.0S  Inches,  the  length  of  each  side 
12.0  Inches  (in  plan-form),  E • 30  x 106  lb/in2, 

Poisson's  ratio  * 0.3  and  the  material  density  is 
— 3 p 4 

0.732  x 10  lb-secVin  . Boundary  conditions  to  be 
considered  are  (a)  all  sides  simply  supported,  and 
(b)  all  sides  clamped.  Natural  frequencies  and 
associated  mode  shapes  are  desired.  The  natural 
frequencies  are  indicated  by  Equation  (37)  and 
a computer  program  for  determination  of  these  fre- 
quencies, associated  mode  shapes,  and  boundary 
residuals  is  available  upon  request  from  the  authors. 

For  the  case  of  simply  supported  edges  the  first 
four  natural  frequencies  are  indicated  below.  Results 
are  compared  to  those  due  to  Vlasov  [5]  as  well  as 
Malkina  [6].  The  investigation  [6]  applies  to  the  case 
of  a spherical  dome  with  circular  plan-form  but  is 
employed  here  as  an  approximation  by  considering  the 
greatest  circular  base  to  be  inscribed  in  the  square 
plan-form  shell  under  consideration. 


Frequency 

Present  method 
L*1 
L-7 

Vlasov  [5] 
Malkina  [6] 


1 

6748.86 

6748.85 

6749.51 

6749.42 


6759.30 

6759.30 

6759.29 


6763.15 

6863.63 


6829.48 

6795.35 

6827.35 


6762.31 


The  computer  program  developed  indicates  the  small 
residual  bending  moments  along  each  of  the  sides  of 
the  shell.  For  example,  for  the  first  natural 
frequency  if  the  peak  bending  moment  at  the  mid-point 
of  the  shell  is  taken  to  be  unity,  the  root-mean- 
square-  boundary  residual  bending  moment  (for  L*7)  is 
found  to  be  0.37  x 10‘2.  This  same  function  for  the 
fourth  natural  frequency  is  found  to  be  0.88  x 10"5. 
Clearly  these  constitute  very  adequate  satisfaction  of 
boundary  conditions. 


Again,  the  computer  program  displays  the  small  residual 
deflections  along  the  shell  boundary.  For  example,  for 
the  first  natural  frequency  the  root-mean-square 
boundary  deflection  is  0.99  x 10"3  (for  L=3 ) and  for 
the  fourth  natural  frequency  it  is  0.75  x 10~5  (for 
L=3)  compared  to  peak  mid-point  deflection  taken  to  be 
unity.  Mode  shapes  are  found  by  calculating  the 
relative  deflections  of  closely  spaced  points  on  the 
shell  then  connecting  all  zero-deflection  points  at  a 
given  frequency. 


DISCUSSION  AND  CONCLUSIONS 


The  present  investigation  has  indicated  that,  for 
the  geometry  considered,  the  Edge  Function  Method  is 
well-suited  to  determination  of  natural  frequencies 
and  associated  mode  shapes  of  thin  elastic  shallow 
shells  undergoing  free  vibration.  For  the  cases 
discussed  through  specific  examples,  the  present 
analysis  yielded  results  in  excellent  agreement  with 
existing  analytical  results. 

The  computer  program  developed  durinn  the  course  of 
the  present  investigation  is  applicable  to  any  spherical 
shell  of  n-sided  plan-form.  Further,  the  boundary 
conditions  along  each  edge  are  completely  arbitrary  and 
may  be  different  on  the  various  edges  without  giving 
rise  to  complications  in  frequency  determination. 

The  use  of  the  present  technique  of  analysis  for 
determination  of  natural  frequencies  and  mode  shapes 
of  thin  elastic  shells  yields  an  additional  bonus  not 
orginarily  found  in  other  approximate  methods,  namely 
rapid  determine tioli  of  how  well  the  specified  boundary 
conditions  have  been  satisfied.  This  is  accomplished 
through  digital  evaluation  of  significant  structural 
parameters  along  the  boundary,  say  deflection  and  slope 
for  a clamped  edge,  which  should,  of  course  vanish 
identically  but  the  small  residuals  corresponding  to 
each  of  these  parameters  are  routinely  printed  out  by 
the  computer  and  root-mean-square  values  along  any  of 
the  polygonal  edges  displayed.  Thus,  one  may  readily 
ascertain  exactly  how  well  the  specified  boundary 
conditions  have  been  satisfied  along  all  boundaries  of 
the  plate. 
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APPENDIX 

SHELL  POLAR  FUNCTIONS 
Equation  (9)  gives  for 

V V - p)0  > 0 (A-l ) 


one  solution  of  which  is 
a 

V 0 * 0 

from  which,  following  O'Callaghan  [4] 

4 * Re[Ezx  + FzzX+1];  z = x + 1y  (A-2) 

where  E and  F are  complex  constants  and  X is  arbitrary. 
Accordingly  on  substituting  for  V40in  (7),  we  obtain 


V2W  ■ 0 

and  equation  (6)  then  gives 

W * V2* 

phi/ 

On  introducing  the  operators 

_3 
3x 


10  3 -10  3 

e R 57  + e g- 


3 , ,J0„  3 < 14>  3 

V 31  qi? 


(A-3) 

(A-4) 

(A-5) 


from  which 

32  32  32 

7 7 iff 

and  on  operating  with  V2  on  (A-2)  and  substituting  in 
(A-4)  we  obtain 

« ’ Re[FzX]  * Re(F'zX]  (A-6) 


where 

F*  * qF;  q . (A-7) 

ohru 

To  fit  into  the  general  polar  form  for  the  tth  derived 
function 

Qt  - Re[EAzX*  ♦ F(BzzX*+1  ♦ CzX*+2)]  (A-8) 

where  A,  8 and  C are  complex  functions  of  0 and  X,  and 
where  necessary  will  have  a subscript  t - A{,  B{,  Ct  - 
to  distinguish  the  different  functions.  Equation  (A-6) 
for  W then  gives: 


A ■ 0;  8 • Oi  C * 


Values  of  A,  B,  and  C for 
parameters  0,  u^,  vq,  etc 


function  t 

X*  At 

W 

1 

X-2  q 

0 

2 

X 1 

Uq 

3 

X-l  A0e^q 

vq 

4 

x-i  ia3 

®q 

Nq 

5 

6 

X-l  A4/r 
X-2  (X-l)e21*i 

Mq 

7 

X-4  0 

Mq 

"nt 

8 

X-4  0 

Tq 

9 

X-2  A12 

Vq 

10 

X-5  0 

11 

X-l  A 4/r 

■fit 

12 

X-2  -1X(X-l)e‘ 

!;  X » X-2  (A-9! 

each  of  the  remaining  eleven 
are  tabulated  below. 

Bt  ct 

0 0 

1 0 

■o*1* 

iB3  -1Cj 

B4/r  C4/r-1  C^fyz2 

(X+l )e21*q  2(X+1) 

0 D*X(X-1 )e21*q 

0 1C; 

B12  clZ+cp/rz2 

0 (X-2)i*q  Cg/r 

B4/r  yr-Xe^q  C,/z2 

*q  -1XX'  0 


0*  * D(l-v) 

X'  - (X+l)e21*q 


' 


